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ABSTRACT 

Until recently, IR observations of comets have consisted mainly of low 

resolution photometric measurements of grain (dust) radiation. The present 

paper explores the possibility of detecting IR molecular line emission from 

cometary parent molecules. Due to the non-LTE conditions in the inner coma and 

the large amount of near IR solar flux, IR fluorescence will be a significant 

source of cometary emission and, in fact, will dominate the grain radiation in a 

sufficiently high resolution instrument. The detection of this line emission 

will be difficult due to absorption in the terrestrial atmosphere, but it 

appears possible to measure cometary 1^0 emission from airplane altitudes (^12 

km). As IR molecular line emission represents one of the few promising methods 

» 

of detecting cometary parent molecules directly, further research on this 
problem should be vigorously pursued. 

Subject headings: comets — infrared: spectra -- molecular processes 


PRECEDING EAjGE jBL'ANK NOT FEMEEf 


I. Introduction 
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Infrared (IR) observations of comets have generally been made using filter 

«* 

photometry. Acquired at low spectral resolution (x/ax £ 100), these studies 
provide data primarily on the cometary grains (dust), which are small (typically 
micron size but sometimes larger) solid particles dragged away from the nucleus 
by the subliming "parent" molecules (constituents of the nucleus which leave in 
the gas phase). Except for the so-called "silicate features" at *10 ^m and *18 
^m, the IR radiation from the grains Is primarily thermal (i.e., greybody In 
nature), so that photometry Is well suited for the study of grain properties. 

IR photometry of bright comets is reviewed by Ney (1982), and of faint periodic 
comets by Campins and Hanner (1982). 

This paper is concerned with IR cometary radiation which is of a completely 
different nature from that described above, and whose detection and 
characterization will, in general, require observations at a higher spectral 
resolution than filter photometry can provide. We investigate the possible 
existence of IR molecular line emissions in comets. In particular, we 
^investigate whether or not there will be detectable IR molecular emissions from 
the cometary parent molecules. Since many of the proposed parent molecules have 
strong ro-vibrational transitions in the IR (e.g., 1^0, C0 2 , CO, NH^, CH^, • . • . ) * 
and since some of these do not fluoresce in the visible or ultraviolet (UV) 
spectral ranges (e.g., H 2 0, C0 2 , NH 3 , CH 4 ,...)» and others do not have electric 
dipole allowed rotational transitions (e.g., CO^, CH^,...), the IR becomes a 
prime candidate spectral region for the discovery and study of these molecules. 

<i * 

IR molecular emissions have already been detected in comet West (see 
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Johnson, Fink, and Larson 1983), but these emissions were from electronic 
transitions or' "fragment" molecules (molecules which do not reside, as such, in 

;v 

the nucleus but are created in the coma, for example, by photodissociation of 
parent molecules)., The physics of these emissions is therefore very similar to 
the physics of visible or UV cometary emissions. In contrast, we discuss the 
production of IR molecular emissions due to ro-vibrational transitions in 
cometary parent molecules. The molecules we have chosen to investigate are H^O, 
COg, CO, and NH-j. These molecules do not form a complete set of possible 
cometary parent molecules, but they serve the purpose of illustrating basic 
principles presented in this paper. Furthermore, since is generally 
considered to be the most important (i.e., most abundant) molecular constituent 
of the cometary nucleus, while CO^ and CO are regarded as particularly important 
in some "new" (i.e., long period or non-periodic) comets (see Del semme 1982 for 
a review), we have attempted to be realistic in our choice. Although NH 3 is no 
longer considered to be very abundant in the nucleus (<1% by number), it was 
proposed as an important constituent in early comet models (Whipple 1950 and 
1951) and continues to be used as a constituent of hypothetical models of 
cometary nuclei (Huebner and Giguere 1980; Mitchell, Prasad, and Huntress 1981; 
Swift and Mitchell 1981; Biermann, Giguere, and Huebner 1982). Laser heterodyne 
spectroscopy may be a sensitive technique for the detection of cometary NH 3 , and 
since we will use this technique to search for NH^ in bright comets, we have 
included' NH 3 in the list of molecules to be considered here. 

We have also decided to include ro-vibrational transitions of OH in our 
study. Although OH itself has never been proposed as a constituent of the 
nucleus, its presence in the coma is generally considered to be a tracer for 
H^O. OH has already been observed in many comets in the UV and radio regions 
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" (see Feldman 1982 and Despois et al . 1981 fop a discussion), so It is 

Interesting to consider its possible IR emissions. Numerous fluorescence 

calculations have been made for OH to interpret Its UV and radio emissions 

(Hunaerts 1953; Dossin et al . 1961; Blraud et al. 1974; Mies 1974a; Lane, 

Stockton , and Mies 1974; Despois et al . 1981; Schleicher and A' Hearn 1982). The 

« 

most recent calculation (Schleicher and A' Hearn 1982) appears to be the most 
accurate since it is based upon a significantly improved solar spectrum and more 
molecular transitions, compared with earlier formulations. However, none of 
these previous works have included IR pumping and thus they may significantly 
underestimate the amount of IR emission. Part of the motivation for including 
OH in the present study is to demonstrate the importance of IR pumping in 
producing excitation in cometary molecules. 


This paper is organized as follows. First, we consider in a general way 
the physical processes operating in the inner coma of a comet. Parent molecules 
are most abundant in this region and, hence, the inner coma will be the primary 
source of their emission. Next, we discuss specific molecules to determine the 
regions in the coma where local thermodynamic equilibrium (LTE) is not valid. 

We then demonstrate that scattering of solar radiation can be a significant 
source of [R molecular line emission when LTE is not maintained. For this case, 

we discuss excitation rates for IR fluorescence and estimate the brightnesses of 

\ 

IR molecular line emissions from a representati ve moderately-bright comet. 
Finally, we summarize our results, discuss their significance, and comment on 
the future role of IR observations of comets. 

Concerning previous work in this area, a literature search reveals that a 
virtual void previously existed. The present paper is actually an extension of 
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earlier work by one of us (Mumma 1980), which to our knowledge Is the first 
attempt to explain IR cometary molecular emission by addressing the non-LTE 
effects which are principally responsible for these emissions. Recently, 

i. 

Yumamoto (1982) presented calculations for the fluxes expected from various 
possible cometary molecules in the IR. His paper is a useful treatment of this 
subject, but he gives virtually no details concerning his calculations and he 
generally neglects many of the points considered in the present paper. For 
example, Yamamoto assumes that molecular line radiation in the coma is due to IR 
fluorescence, implying that vibrational LTE is not valid. On the other hand, he 
assumes that the lines within the vibrational bands follow a thermal 
distribution, implying that rotational LTE is maintained. Since neither 
assumption is obviously valid throughout the coma, each requires detailed 
discussion. Even more recently, a French group (Encrenaz et al . 1982; Crovisier 
1982) has presented, in a preliminary fashion, some results of an analysis they 
have performed to estimate the IR emissions from comet Halley. The only 
molecule for which they give any detail is HgQ, and, in this case, their results 
appear to be consistent with our own. 

II. General Considerations 

Fig. 1 shows the IR flux from comet Bennett (1970 II) measured by Ney 
(1974) using filter photometry. As discussed earlier, the dominant component of 
this observed flux is thermal radiation from cometary grains. Also indicated in 
Fig. l are the positions of the vibrational transitions considered in this 
paper. The presence of a strong internal source of IR radiation raises the 
interesting possibility of detecting parent molecules in the IR by their 
absorption of thermal dust continuum. However, collisional relaxation does not 
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compete favorably with radiative vibrational relaxation in the cometary coma, 
Implying that each photon absorbed from the thermal dust continuum is 
re-emitted (see Sec. III). Since both grains and gas are Isotropic emitters, 
and since they have similar spatial distributions, the surface brightness is 
unchanged In the spectral line. Neglecting branching ratio effects, it appears 
that cometary parent molecules cannot be detected by absorption lines In the 1R. 

Cometary atoms and molecules are detected in visible and UV spectra by 
their emission lines, which are usually produced by fluorescence in the solar 
radiation field ([01] is an important exception). The solar flux (Fig. 1; Labs 
and Neckel 1968) actually peaks in the near IR (\ ^1.6 um) suggesting that IR 
fluorescence may provide a means of detecting cometary parent molecules if they 
re-emit the absorbed solar radiation before collisional (or other) de-excitation 
occurs. However, it is not obvious that fluorescence will be important for IR 
transitions since the competition between radiative and collisional processes is 
much keener there than in the UV or visible. Collision cross sections generally 
grow for decreasing energy separations of the levels involved in transitions, 
while spontaneous radiative transition probabilities usually decrease under the 
same circumstances. Thus, even though collisions are unimportant tor most UV 
and visible transitions, they may still play an important role in some IR 
transitions, and this role must be determined individually for each molecule. 

If LTE is maintained in the coma, then IR fluorescence in the solar 
radiation field will not occur, because of quenching by collisions. Under these 
circumstances, the IR radiation from the molecules in the coma will be primarily 
thermal in nature. This is not necessarily a drawback with regard to the 
observation of IR molecular line emission. For example, a simple calculation 
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* demonstrates that the v 2 level of C0 2 (which Is 667 cm" 1 above the ground level) 
would have a significantly higher population under the conditions of vibrational 
LTE at a temperature of 200 K than it would have if the level were in 
fluorescence equilibrium. The reverse situation prevails (i.e., larger 
population In fluorescence equilibrium) for levels having larger excitation 
energies (which prevents easy thermal excitation) and for which a larger amount 
of solar flux is available for fluorescence. In fact, for all of the other IR 
bands considered in this paper, the IR emission will be significantly enhanced 
when vibrational LTE is not maintained. 

The relative importance of colli si onal and radiative effects may be 
assessed by examining the rate equation for an excited energy level. A general 
rate equation taking into account all possible physical processes would be 
hopelessly complicated. Fortunately, considerable simplification is made 
possible by eliminating those effects which are expected to be unimportant in 
the cometary coma. 

i. 

» 

We will not consider a time-dependent rate equation. Following Arpigny 
(1964), the excitation and de-excitation rates considered here will generally be 
large enough that the molecules will at least approach the steady state 
condition. A probable exception is CO, due to its relatively long rotational 
radiative lifetimes. Preliminary results for CO indicate that the relative 
vibrational populations are equilibrated very quickly, but the distribution 
among the rotational levels may never reach steady state in the observable coma 
(Weaver and Chin, private communication (1983)). The rotational distributions 
used for our calculations are discussed in Section IV. 
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Excepting homonuclear molecules (e.g., C 2 , H 2 ), we expect nearly all 
molecules to be in the lowest vibrational level of the ground electronic state. 
Thus, for both colli si onal and radiative processes, we assume that at least one 

r 

of the interacting levels is In this state. 

Restricting interactions to those connecting the excited level and the 
ground state only, eliminates from the discussion a mechanism we term 
"Franck -Condon pumping," which can produce significant IR emissions. In this 
process, molecules excited by resonance electronic transitions in the visible or 
UV spectral regions may be left in higher vibrational states upon returning to 
the ground electronic state. These high vibrational levels will subsequently 
decay by vibrational cascade, accompanied by the emission of IR radiation. For 
favorable cases many IR photons will be emitted for each UV or visible photon 
absorbed, thereby producing larger photon emission rates in the IR than in the 
UV or visible. In general, a mismatch in the mean bond-lengths for the upper 
and lower electronic states will result in a Franck-Condon envelope which is 
extended over several vibrational levels, and vibrational pumping will result 
for such molecules. For example, the CN (A ^ir - X red system has a wide 
Franck-Condon envelope and fluorescence in this system will produce extensive 
vibrational pumping, and, consequently, large IR emissions. Although 
Franck-Condon pumping works well for most diatomic molecules, it will not be as 
important for polyatomic molecules because their upper electronic states usually 
pre-di ssociate. 

Stimulated emission and optical trapping effects are neglected in the rate 
equation. A quick calculation demonstrates that stimulated emission is 
certainly unimportant for the vibrational transitions considered, and is 
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probably unimportant for the rotational transitions as well. In the inner coma, 
where the column densities are highest, fairly large optical depths (>>l) may be 
encountered in the cores of the strongest ro-vibrational transitions of some 
parent molecules. When significant radiative trapping occurs, a resonance 
photon may be scattered many tin$a$ before exiting the coma, and it has a 
correspondingly higher probability of being thermal i zed. This phenomenon lowers 
the "effective" Einstein A value for the transition and allows LTE to be 
maintained at a lower density than would have been required in the optically 
thin case. This effect plays an important role in the analysis of non-LTE 
processes in the mesospheres of Mars and Venus (Deming and Mumma 1983; Dickinson 
1972, 1976) and may be important in comets. We include optical depth effects 
later in the discussion of IR fluorescence (section IV), but we will neglect 
optical trapping in determining where in the coma the assumption of LTE is not 
v a 1 id. 

Only neutral -neutral collisions producing energy exchange between 
translational and vibrational (T«/ process) or rotational (T-R process) modes 
are considered. Collisions transferring vibrational energy of one species into 
vibrational energy of the same or other species (V-V process) may be important 
for high excitation energy vibrational modes, but we will not consider this 
process. Excitation and de-excitation of neutral molecules can also occur by 
collisions with ions and electrons (electrons being the most important-). 

Although fairly high electron densities may be present in the inner coma (n e £ 
10^ cm"** according to Giguere and Huebner 1978), we feel that not enough 
quantitative information is available to consider this process realistically 


here. 
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We also neglect direct production of molecules in excited states. For 
example, some parent molecules leave the nucleus with thermal excitation and may 
relax radiatively as they flow Into the coma. Calculations indicate that this 
will not be a significant source of emission for the transitions considered 
here. It is also possible that CO is a photodissoclatior. product and, thus, may 
be produced in an excited state (e.g., in the photodissociation of C0 2 (Huebner 
and Carpenter 1979)). Highly excited rotational states of OH are produced in 
the photodissociation of H 2 0 by H Ly a photons (Yamashita 1975). Since we are 
concentrating on parent molecules, we neglect such effects. 

Thus, we write the following simplified rate equation: 

d V dt - 0 ■ n t (B lu »+C t(J ) - n u (A^+C^) (1) 

where, n , n M s number densities in the lower and upper energy levels, 
i u 

respectively, of the effective two level system, 

p 3 density of the radiation field at the position of the 
molecule and at the frequency of the transition, 

B = Einstein coefficient for stimulated absorption, 
tu 

A 1( 2 Einstein coefficient for spontaneous emission, 

ug, 

C . C = collisional excitation and relaxation rates, respectively, 
tu* ug, - 


and 
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Equation (1) is an acceptably accurate and tractable description of the f 
physics in the inner coma, enabling us to discuss our ideas concerning the 
production of IR molecular emissions quantitatively. Processes not included 
explicitly are assumed to produce relatively small effects and can be Included 
in future, more detailed studies. 


We next consider (separately) vibrational transitions in the near to middle 

IR and rotational transitions which occur at much longer wavelengths. We treat 

the former in order to estimate the integrated IR brightnesses of vibrational 

bands. We treat the latter because these transitions help determine the 

population distribution among the various rotational levels of the ground state, 

and, hence, the intensity distribution within the vibrational bands. For 

vibrational transitions, n will represent the lowest energy vibrational level, 

« 

while n u will be an excited vibrational level. For rotational transitions, n 4 
and n u will refer to energy levels within the lowest energy vibrational level 
(for radiative transitions this will mean energy levels with total angular 
momentum quantum numbers differing by one unit). 




The collisional excitation and relaxation rates are given by: 


C *u 5 n T° t v 


( 2 ) 


and C Ul = n T o f v 


( 3 ) 


where n^ = total gas number density (cm’ 3 ). 


0 t O 

t > 


excitation and relaxation cross sections (cm 3 ), 
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respectively, for the (T-V) or (T-R) collisions 

connecting the levels considered, 

and * i average relative speed of the colliding ©alltcules. 

Using the detailed balance condition, equation (1) can be rewritten as; 

V\ * (9 u /9 l )(e-( hc ' ,/ ' tT) )C(l + B tU p/C w )/(l+A Ul /C Ut )3 ' (*) 

where, T 3 local kinetic temperature (K) , 

and v a wavenumber of the transition (cm"*). 

From equation (4) we see that a sufficient condition for LTE is; 

B iu p/C iu «1 and A ut /C ul «1 

In other words, LTE is maintained in the inner coma when collisions dominate 
radiative processes (Notice that setting U +b * u p/C 4 J*(lM Ui /C ut ) yields another 
sufficient condition for LTE. However, this condition simply implies that LTE 
is valid when the radiation field is identical to that of a blackbody whose 
temperature is equal to the kinetic temperature of the coma). 

It is clear that radiative excitation plays an important role in the 
excitation of vibrational transitions. However, the different spectral 
character of the solar and grain radiation fields raises the interesting 
question: which source dominates the local radiation field in the 10 urn 
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spectral region at the position or' the comet? The temperature of a typical 

r 

cometary grain at R (s heliocentric distance of comet) * 1 a.u. Is *330 K 
(following Ney (1974), this Is *20% larger than the temperature of a blackbody 
at this distance fra the sun). The radiation density near 10 ^m of a blackbody 
at T ■ 330 K Is approximately three orders of magnitude larger than the 
corresponding density of the solar radiation field at. 1 a.u. However, the 
radiation density due to the grains will not be equal to that of a blackbody at 
the same temperature as the grains for two reasons. First, the grains are not 
"black"; they have an emissivity less than one and their emissivity changes with 
wavelength (Ney 1974). More importantly, since the grains do not completely 
fill the region around the nucleus (i.e., they form a very tenuous 
"atmosphere"), the grain radiation field is diluted from the blackbody value. 

For example, Ney (1974) estimates for comet Bennett (1970 II) that the average 
grain optical depth is *10“ 4 -10~ 5 in a field of view subtending *7000 km at the 
comet, producing a dilution factor of this same order. We calculate that for 
wavelengths longward of *lOum, the grain radiation field will be stronger than 
the solar field for a region within *500 km of the nucleus. Nevertheless, for 
the IR vibrational transitions considered here, we consider only fluorescence 
produced by the solar field. Under these conditions we have: 


V * %and 


where $band 2 excitation factor or "g-factor" for the vibrational 

band (photons s’ 1 molecule" 1 ). 

The concept of the g-factor is extensively discussed in section IV when we 
consider IR fluorescence. 
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If the total number density is taken to be that of fUO and a Haser model 
(Haser 1966; the validity of the Hassr model is discussed by Festou 1981) is 
used to determine the HgO density, then: 

n T - (^H 2 Q/ 4irr2v ') e ^ r/V ' T ^ (S) 

where, q H Q production rate (molecules s’ 1 ), 

r = radial distance to the center of the nucleus, 

v' * radial outflow velocity tor the HgO molecules (assumed to be 
the same value for all HgO molecules), 

and t = H 2 0 lifetime. 

For H,pO , v't *4-8 X 10 4 km at R = 1 a.u. (Weaver et al . 1981) so that the 
exponential term may be safely neglected within a few thousand kilometers of 
the nucleus. 

We now define the "disequilibrium distance", 0 as the distance from the 

r 

nucleus outside of which collisions will be of minor importance in producing 
vibrational excitation. We calculate by setting the radiative excitation 
rate equal to the collisional rate. Using equations (2). (5), and (6) (but 
neglecting the exponential factor in (6)), the disequilibrium distance can be 
written as: 


( 7 ) 
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The temperature of an H^O-lce nucleus changes very little with heliocentric 
distance when R < 2 a.u. because virtually all the incident solar energy goes 
into sublimation and little is available for surface heating (Weissman and 
Kieft'er 1981). The temperature attained by the nucleus depends on factors such 
as albedo, heat transport parameters, etc., but "reasonable" values for all 
these parameters yield temperatures of *200 K. The kinetic temperature of the 
sublimed gases will initially decrease with increasing distance from the nucleus 
due to expansion. Photodissociation heats the coma, but probably not by enough 
to overcome the expansion cooling (at least in the inner coma). These 
conclusions are supported by the calculations of Shimizu (1976) and Marconi and 
Mendis (1982). Besides expansion cooling and photochemical heating, these 
authors have also included radiative cooling by pure rotational transitions of 
H 2 0. They find that the kinetic temperature decreases from a value of *200 K at 
the surface of the nucleus tc *10 K at *100 km from the nucleus, followed by a 
monotonical ly rising temperature for larger distances from the nucleus. 

However, we will assume that the temperature of the nucleus (T a 200 K) also 
defines the kinetic temperature of the gas throughout the inner coma. We simply 
note that the processes not considered here will probably shrink the 
disequilibrium distance to a smaller value than we calculate. Furthermore, by 
analogy with work done for interstellar clouds (Nuth and Donn .1981) , 
considerable departures, from LTE can be expected well within the radius at which 
the colli si onal and radiative rates are equal. 


Under our assumptions, the average relative thermal speed of the colliding 


coma molecules is given by: 
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where, k s Boltzmann's constant, 


and u 3 reduced mass of the colliding particles. 


Since we assume that HgO is the most abundant molecule, the coma can be 
considered as a mixture of molecules In an H 2 0 gas, i.e., H 2 0 is always at least 
one of the collision partners. The values of the average relative speeds for 
HgO collisions with H 2 0, CU 2> CO, and at T*200 K. lie between 0.57 km s" 1 and 
0.73 km s' 1 . All of these are fairly close to the average outflow velocity of 
the gas (Delsemme 1982), so we will assume v * v‘. Then, equation (7) 
simpl ifies to: 

° f a ^H 2 0° + /4irg band) 1/2 

In a similar manner a disequilibrium distance, D + , can be defined as the 
distance from the nucleus outside of which vibrational radiative relaxation 
begins to dominate collisional de-excitation. Setting the radiative relaxation 
rate equal to the collisional relaxation rate yields: 

- <y H2 oV*’V 1/2 


( 10 ) 
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' It is important to discuss disequilibrium distances for both excitation and 

« 

relaxation because it is possible, for example, for collisions to dominate the 
relaxation of the molecule while being unimportant in its excitation. If this 
occurs then collisions will quench the potential IR fluorescence. 


For rotational transitions, the radiative pumping rates should be very 
small because both the radiation densities and the transition probabilities are 
much smaller than their values for electronic and allowed vibrational 
transitions. In this case, we see from equation (4) that LTE holds only when 
A << C . Setting the colli si onal and radiative relaxation rates equal to 
each other defines an expression for the rotational disequilibrium distance. 

This expression can be obtained from equation (10) above by substituting 
quantities appropriate for pure rotational transitions for the corresponding 
vibrational quantities. 

When colli si onal effects are unimportant, allowing rotational relaxation, 
the number of lines in a particular vibrational emission band will be greatly 
reduced but the brightnesses of the remaining lines will be greatly enhanced. 
However, fluorescence in vibrational or electronic transitions may result in 
significant rotational pumping, reducing this enhancement. In the extreme case, 
the effective rotational temperature could be as high as the color temperature 
of the sun (this effect has been observed, for example, in C 2 ), Even when we 
find that the rotational disequilibrium distance derived using equation (10) is 
small, we must still assess whether complete rotational relaxation does indeed 
occur by comparing the Einstein A's for rotational transitions to the g, factors 
for vibrational and electronic transitions. 


'III. Individual Molecules 
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The disequilibrium distance Is not very sensitive to the specific value 
chosen for 0 , as long as a "reasonable" number Is used, i.e. , the number 
should be representative of observed comets. Table 1 lists water production 
rates at R « 1 a.u. for several comets which range from "moderately-brlght" to 
"bright". We choose ^ ■ 2X10 29 molecules s’ 1 as an appropriate value. 


Collisions having HgQ as a partner are unusually efficient, and even (T-V) 
cross sections can approach the gas kinetic values. Table 2 lists the collision 
relaxation cross sections ( 0y+ ) for (T-V) collisions between HgU and HgO, CO,, 
NH 3 , and CO. These values were derived using experimental data from Taylor and 
Bltterman (1969) and Hi rsehf elder , Curtiss, and Bird (1964) scaled to the 
temperature of interest (T * 200 K). Since the experimental data usually 
covered a fairly wide range of temperatures and extended down to T^300 K, we 
feel that our extrapolation to 200 K could not be significantly incorrect. For 
NH 3 , we have made the additional assumption that NHg-HgO’col 1 isions are similar 
to NHg-NHg collisions. Also, for CO we assume that Ng-HgU or collisions are 
indicative of CO-HgO collisions (according to Taylor and Bitterman 1969, Mg and 
CO are expected to behave similarly). The cross sections given refer to the 
lowest energy vibrational fundamental for each molecule. Collisions relaxing 
more highly excited vibrational levels must exchange more energy and will 
generally have smaller cross sections than those listed (However, see the 
discussion in section II concerning the role of (V-V) collision processes). 


Rotational transitions (T-R) are expected to occur for nearly every 
collision. Experimental data on collisional rotational relaxation are difficult 
to find and are probably non-existent for the specific interactions considered 
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here. Generally, those data which do exist indicate cross sections comparable 
to the gas kinetic values, so we have assumed that the (T-R) cross sections 
(o^ R ) are identical with the gas kinetic values. These have been calculated 
using data from Hirschfelder, Curtiss, annd Bird (1964) and are listed in Table 
2 . 

t 

Radiative pumping rates are required for calculating the disequilibrium 
distance (D + ) for vibrational excitation (eqtn. 9). Chamberlain ( 1978) defines 
the "g-factor" to be the fluorescence rate due to scattering of the local solar 
flux. For a two-level system, 

g * (*e 2 /mc 2 ) f F Q (11) 

where *e 2 /mc 2 * 8.82 X IQ' 13 cm, 

f = absorption oscillator strength for the transition, 
and F e * the solar flux at the position of the comet, which is 

assumed to have R » l a.u. 

(photons cm " 2 s'Vcm’ 1 ) 

In IR spectroscopy, the “absorption strength" (S line , cm molecule' 1 ) is commonly 
used instead of the oscillator strength (see Penner 1959 or Pugh and Rao 1976). 
The two quantities are related by: 

. % 

s 1ine ' (%e 2 /">c 2 W Mne q (12) 

q = fraction of the total number of molecules which are in the 


where , 
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absorbing level. It is computed in the usual manner from 
Boltzmann factors and the partition function as described, for 
example, in Herzberg (1945). 


We generalize the above expressions to a vibrational band by writing (the solar 
flux does not vary appreciably within each vibrational band): 


9 band * ^0^1 ine * ^Q^band 

The band strengths for H 2 Q and CO are taken from Rothman (1981), and those for 
COjj are taken from Rothman and Benedict (1978). The band strength for the NH 3 
V £ transition is taken from Taylor (1973). Solar fluxes are obtained from Labs 
and Neckel (1968). The band absorption strengths, band g-factors, and the solar 
fluxes used are given in table 3 (see also Sec. IV). 

Using the collisional cross sections (Table 2) and the band g-factors 
(Table 3), we calculate (eqtn. 9) vibrational disequilibrium distances for 
excitation (0 , Table 4). Radiative excitation will apparently dominate 

t 

collisional excitation throughout most of the cometary coma. The largest 
disequilibrium distance is for the \> 2 bancl oT CU 2’ but even in Case 
scattering of solar photons will be the most important excitation mechanism for 
regions farther than *400 km from the nucleus. 

The radiative transition rate (A line , s -1 ) for a ro-vibrational transition 


is given by: 


LC. 


A 11ne ’ (Swcd/q) v 2 S ]1ne 
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( 14 ) 


where, c s speed of light in vacuum (cm s' 1 ). 


and d * ratio of the statistical weight of the lower state 
to that of the upper state. 

Thus, we define the band Einstein A value (A band , s’ 1 ) by: 

2 

A band ‘ ® TCdv band s band < 15 > 

where, "band 5 wavenumber r'er the band origin (cm* 1 ). 

I ' 

The Einstein A's and wavenumbers for all of the vibrational bands considered in 
this paper are given in Table 4. 

Using the Einstein A's (Table 4) and col 1 i si ona 1 cross sections (Table 2), 
we calculate (eqtn. 10) vibrational disequilibrium distances for relaxation (D , 
Table 3). Despite the relatively large values tor some of the colli si onal cross 
sections, all the D's are very small. Apparently, vibrational disequilibrium 
will occur very close to the surface of the nucleus for all the molecules 
considered. 

Rotational disequilibrium distances for relaxation are calculated using 
rotational Einstein A's and the cross sections for (T-R) collisions from Table 
2. Due to the dependence of the Einstein A's, their values can vary over a 
few orders of magnitude depending upon whether we are referring to low J or high 
J (J = total angular momentum quantum number) transitions. The A's listed in 
Table 3 are typical low J values (values for higher J transitions are larger). 
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I 

- We calculate the Einstein A's tor NH 3 and CO following Somerville (1977), and we 

I 

use the values given for HgU by be Jong (1973). The derived rotational 
disequilibrium distances are given in Table 4, and their significance is 
discussed in Section IV. 

Disequilibrium distances are not calculated for OH. Since OH is a jj 

I 

dissociation product, it is round mainly in the outer coma, where densities are 
iow. Also, vibrational and rotational transitions are electric dipole allowed 
in OH so that it relaxes rapidly after excitation. For these reasons, 
collisions are expected to be relatively unimportant for this molecule. This 
conclusion is apparently supported by the agreement between the theoretical 
results, derived assuming pure fluorescence equilibrium in OH, and the 
observations of OH in comets (Schleicher and A'Hearn 1982). 

IV. IR Fluorescence 

1 

* 

Having shown that collisions cannot maintain vib ational LTE for the ) 

1 

molecules considered, we can now calculate their IR fluorescence rates. This \ 

S 

l 

process will be the primary source of IR photons from the cometary molecules. f 

! 

1 

I 

I 

The line-integrated volume emission rate (j) can be expressed as: j 

I 

j =* g n(r) (photons cm" 3 s" 1 ) (16) 

where, n(r) = the number density of the molecule being considered at the -- 

radial distance "r" from the nucleus (cm" 3 ). 
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* Under optically thin conditions, the surface (or column) brightness (B) is given 
by the integral of the volume emission rate along the line of sight (Chamberlain 
1978). Assuming that 'V does not vary with position in the coma; 

B ■ 4* I * g N ( 17 ) 

where, B * surface brightness (photons cm" 2 s"*, usually 

expressed in Rayleighs where; 1 Rayleigh ■ 10 6 photons 
cm“ 2 s“ l ), 

L » line Integrated specific intensity (photons cm" 2 
s" 1 sr' 1 ), 

and N = column density of the molecule in question (cm“ c ). 

Our assumption that "g" does not vary along the line of sight is equivalent 
to assuming that each volume .ement in the coma "sees" the same value of solar 
flux. This means that there is negligible optical depth along the path 
connecting the sun with the volume element in question. However, for regions 
very close to the nucleus .large optical depths can be encountered in molecular 
lines for both the path connecting the given volume element to the sun and for 
the path connecting this volume element to the observation point. The radiative 
transfer for the IR fluorescence is thus extremely complicated (Note that proper 
inclusion of the line shape will further complicate matters because no single 
velocity distribution is valid for all points in the coma and the local velocity 
distributions will be neither isotropic nor Maxwellian). The brightnesses 
calculated using equation (17) will be upper limits. Calculating g-r'actors as 
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If optical depth effects could be neglected Is Instructive for several reasons. 
First, when the optical depths are not too large the optically thin g-f actor may 
produce reasonably accurate brightness predictions. Second, outflow effects 
tend to reduce self shielding of the inner coma. Finally, it Is interesting to 
compare the IR g-factors derived here for the cometary parent molecules, to the 
UV and visible g-factors derived for various dissociation products already 
observed In comets. 

To account for optical depth effects in an approximate way, we also 
calculate some IR brightnesses using the following model. Consider a sphere 
centered on the nucleus, whose radius (r ) is chosen such that the line center 
optical depth (t ) measured along the path from the sun to the nucleus, reaches 
unity at the surface of the sphere (the column density at which t Q =i is 
calculated assuming a Doppler line profile characterized by a Kinetic 
temperature of 200 K). We assume that the sphere so defined uniformly scatters 
into 4ir steradians all of the solar flux incident upon it at the center 
frequency of the molecular line. This scattered flux provides the emergent 
intensity at line center for the cometary emission, which is assumed to have a 
rectangular profile with a width of 1.6 km s" 1 . Integrating over the lineshape 
then gives the amount of radiation scattered by the entire coma in the line 
being considered. This procedure neglects details of the fluorescence process, 
but allows an easy calculation of the scattered flux. The physical 
justification for this model is derived from the fact that scattering of the 
solar flux becomes significant at T(J 1, so that little of the solar flux at 
line-center will penetrate deeper into the coma. A rectangular line profile 
crudely approximates the reduction of self-shielding due to outflow. In all 
cases for which we have made calculations this procedure gives a lower 
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* brightness than the optically thin formalism, Indicating that at least we have 
moved In the right direction (although not necessarily to a more accurate 
prediction). 

The molecular densities needed for making brightness calculations are given 
by eqtn. (6), replacing the various parameters pertaining to HgO with those 

appropriate for the molecule under consideration, photodestruction lifetimes 

► 

are given for Individual molecules in Table 5. A radial outflow velocity of 0.8 
km s" 1 is used for all molecules (the actual values probably lie between 0.5 and 
1 km s' 1 ; see Weaver et al . (1981) and Delsemme (1982)). The column density is 
obtained by numerical integration of the density distribution along the line of 
sight. The parameter (r Q ) is obtained by integrating the optical depth along 
the sun -nucleus line, starting at very large values of r (where the density is 
extremely low due to expansion and lifetime effects) and proceeding to the point 
where the optical depth equals unity, at line center. 

Nearly all observational evidence indicates that i s the most abundant 
constituent of cometary nuclei, and so we consider its g-factors first. In 
particular, we consider its v 2 fundamental bang (001-000) at \* 2>66 iim (3756 
cm" 1 ). Fig. 2 is an energy level diagram showing. the lowest few rotational 
energy levels in both the upper and lower vibrational states of this band. The 
molecular ensemble, and thus the energy diagram, is divided into two 
non-interacting parts, the ortho form of H 2 0 having a statistical weight three 
times that of the para form. Angular momentum quantum numbers and the various 
symmetries of the levels are indicated in Fig. 2 (see Herzberg (1945) for 
details). Typical Einstein A ' s and the electric dipole selection rules 
governing these transitions are given in de Jong (1973). Certain rotational 
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* levels form a "backbone'* (de dong 1973), defined by the lowest energy levels for 
a particular value of J (J = total angular momentum quantum number). Using the 
formalism of sections II and III, we find that cometary H 2 Q molecules not 
originally in the backbone levels will very quickly relax to these levels. 

Furthermore, the backbone levels will quickly relax to the lowest possible value 
of J. The relaxation rates decrease as J decreases but, as can be seen from 
Table 4, the rotational populations are expected tc relax within + 1000 km from 
the nucleus even for the lowest energy levels. Thus, the conditions in the coma 
are seen to differ radically from the usual LTE situation. In LTE, many energy 
levels are populated leading to complex spectra (HgU is an asymmetric rotor) 
with numerous lines in each vibrational band. However, in comets we expect to 
see a large reduction in the number of emission lines and, consequently, a 
tremendous enhancement in the brightnesses of those remaining in each 

I 

vibrational band, since all of the bind strength will be concentrated into those f 
few lines . 

We have calculated g-factors for individual lines in the v ~ band of H~0 

4 Cm 

assuming that 25% of H 2 Q molecules are in the 0 Q0 level of the ground state and 
75% of the molecules are in the 1 Q1 level of the ground state. This is the 
expected population distribution when all of the H 2 0 molecules are completely 
relaxed and assuming that the ortho and para forms of H 2 0 have a relative 
abundance equal to the ratio of their statistical weights. Under these 
conditions, only six transitions are possible in the v 3 band, and these are 
indicated in Figure 2. Oscillator strengths and branching ratios were computed 

| 

from the individual line strengths given in the Air Force Geophysics Laboratory 1 

j 

(AFGL) molecular line atlas, using the formalism of Section HI. Thu individual -*i 

line g-factors are given in Figure 3. The g-factor for the band is given in 
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• Table 3 (see Section III). OEPoSr ^Au T y ■ 

Yamamoto (1982) derives a band value which Is 1.13 times larger than ours. 
Since we have not been able to determine exactly how he performed his 
calculation we cannot account for the difference. However, we consider this 
difference to be small and conclude that the two calculations for the band 
g-factor are in agreement. We disagree with Yamamoto (1982) significantly on 
the predicted shape of the band. Yamamoto (1982) assumes an intensity 
distribution characteristic of H^O in LTE at T * 300 K. Our results indicate 
that the number of H 2 0 lines in the v 3 band may be reduced to four very intense 
lines and two weaker lines. We also note that IR pumping will produce several 
intense rotational lines in the far infrared (see Figure 2). 

Encrenaz et al . (1982) derive a band value which is *1.2 smaller than ours. 
Since we use virtually identical values for the IR absorption strengths and 
solar intensity, the difference apparently arises in the calculations 
themselves. The intensity distribution within the band is virtually identical 
in the two cases, reflecting the fact that both groups have independently used 
similar arguments to justify the population distribution described above. 

A component of the cometary nucleus more volatile than H 2 0 and which is 
present in substantial amounts (*10-15% relative to H 2 0), is apparently 
necessary to explain the vaporization behavior of certain new comets. C0 2 and 
CO have both been proposed as candidates for this volatile component (Delsemme 
1982). 

Due to its symmetry, the C0 2 molecule does not have a permanent electric 
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dipole moment and, thus, its rotational levels cannot relax radiatively. In the 
inner coma, the rotational excitation will be dominated by collisions! processes 
while IR fluorescent pumping will determine the rotational distribution in the 
outer coma. A detailed model Is needed to compute the relative contribution of 
each component in a given field of view (FOV). We assume here that the 
intensity distribution within the C0 2 vibrational bands is given by the LTE 
distribution at T*200 K. 

We have calculated g-factors for individual transitions in the \> 3 band of 
C0 2 (00°l-00°0) centered at \ = 4.26 ^m (2349 cm"*), and the results are 
indicated schematically in Figure 3. Each emission line is pumped by two 
absorptions (P and R), and the pumping rates were calculated after scaling the 
AFGL line strength data from 296 K to 200 K. Branching ratios into the emission 
line being considered were derived using eqtn. (14). The v 3 band of C0 2 is one 
of the strongest IR bands known, indicating that CO^ might be detectable in the 
coma to fairly low abundance levels. The band g- factor is listed in Table 4 and 
is approximately a factor of ten larger than the H 2 Q v 3 band g-factor. Yamamoto 
(1982) gives a g-factor for the C0 2 v 3 band which is virtually identical with 
ours. 

The lower rotational levels of CO are maintained in LTE out to fairly large 
distances from the nucleus (Table 4). Since the Einstein A coefficients for 
rotational transitions scale as J 3 , we would expect the higher energy levels to 
relax much more quickly. However, time constants for CO are slow enough that a 
steady state assumption may be invalid (see discussion, section II). A detailed 
calculation addressing CO excitation in cometary comae is in preparation (Weaver 
and Chin 1983). Pending completion of the exact model, we will assume that the 
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•population distribution of the rotational levels Is given by the LTE value at T 
■ 200 K. The intensity distribution In the g-factors for the (1,0) band of CO 
(Figure 3) should be considered as Illustrative only since the actual intensity 
distribution could depart significantly from this. The g-factors were 
calculated In the same manner as described above for C0 2 , using molecular data 
derived from the AFGL compilation. The center of the CO (1,0) band Is at \ • 
4.67 urn (2143 cm" 1 ). Each emission is pumped by two absorptions (P and R). We 
emphasize that direct IR pumping in the (1,0) band was the only excitation 
mechanism retained. This Is also true for the polyatomic iftolecules considered 
here, but other excitation mechanisms may be particularly Important for diatomic 
molecules (cf. Sec. II). The g-factor derived for the (1,0) band is given in 
Table 3. Notice that the band g-factor for CO is very close to the value 
derived earlier for the band of HgO. Yamamoto's (1982) band value for CO is 
in good agreement with our value, differing by only 10%. 

NH 3 is no longer considered to be an abundant constituent of the nucleus, 
but it has special significance for us because we have an instrument (Mumma et 
al . 1982) which may provide a sensitive measure of its abundance. NH 3 is a 
symmetric top molecule which exhibits "inversion doubling" (Herzberg 1945). 
Figure 4 is an energy level diagram for the ground state of NH< 3 and shows the 
inversion doubling, rotational structure, and the symmetry properties of the 
ground vibrational state. Following section II (see Table 4), transitions 
within each "K-ladder" will quickly relax the molecule so that only the d*K 
levels are populated. The transitions with 0 are not allowed under 
electric dipole selection rules, hence the J»K levels are metastable. Since the 
Einstein A coefficients connecting the inversion doubled levels are very small 
(Aj'2 X 10" 7 s’ 1 ), the anti symmetric and symmetric levels form independent 
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- reservoirs within each J * K level (except for J » K » 0 which possesses only an 
antisymmetric level). We thus use a partition function for ccmetary NH 3 in 
which only the J * K levels of NH 3 are populated, and we take the population 
distribution among these levels to be characteristic of the temperature in the 
inner coma, i.e., T * 200 K. . 

We have calculated g-factors for 122 lines in the v 2 band of NH 3 located 
near 10.5 ^m (950 cm* 1 ) using molecular data from Taylor (1973). The values for 
the strong lines are shown in Figure 5. Absolute strengths were calculated 
using the non-LTE partition function, invoking the sum rule, and normalizing to 
the vibrational band strength (Taylor (1973)). Selection rules dictate that 
only Q and R pumping lines are possible. Except when K = 0, a q absorption 
pumps both q and P lines in emission, while an R absorption pumps R, P, and q 
emissions. For K = 0, only an R absorption is possible and this pumps only R 
and P lines in emission. Individual oscillator strengths were calculated using 
equation (12) and equation (14) was utilized to derive branching ratios. The 
g- factor for the v 2 band is listed in Table 3. In contrast to the case of LTE, 
we see that nearly 10% of the entire band emission is concentrated into a single 
line, sq (3,3)1 Yamamoto (1982) gives a g-factor which is a factor of 3.3 
smaller than our value. We do not know the source of this discrepancy. 

Finally, we have calculated g-factors for lines in the (1,0) band of OH 
centered at x = 2.80 u m (3570 cm* 1 ). We assume that all OH molecules are in the 
lowest rotational and vibrational level of the ground electronic state ( ir 3 / 2 )* 
OH is thought to be a photodissociation fragment of H-jO and is found mainly in 
the outer coma where collisions are not important, allowing complete vibrational, 
| and rotational relaxation to occur. We assume that UV fluorescent pumping does 

i! 

a 
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• not produce significant equilibrium populations in its excited vibrational or 
rotational levels. These assumptions are supported by detailed calculations 
which indicate that even in fluorescence equilibrium over 95% of the OH 
molecules are in the lowest energy level (cf. Schleicher and A' Hearn (1982)). 

Under the conditions described above, there are only two possible 
absorptions in the (1,0) band and these pump five emission lines (here a line is 
defined as a sum over both A-doubled components; see Schleicher and A'Hearn 
(1982) for a description of the ground electronic state of OH). Unlike our 
previous calculations we also include the contribution to the (1,0) band 
emission by pumping in the (2,0) band followed by cascade in the (2,1) band. 

The contribution to the (1,0) band g-factors from this excitation mechanism is 
generally about 10-20% of the contribution from direct pumping in the (1,0) 
band. There are two absorptions possible in the (2,0) band and these pump the 
same five emission lines mentioned above plus three more lines (however, since 
these latter three lines are pumped only by absorption in the (2,0) band their 
g-factors are smaller than the others). Einstein A coefficients for the allowed 
transitions are taken from Mies (1976). Although recent work (Goldman 1982) 
suggests that the A's for the (1,0) band transitions should be reduced by a 
factor of 0.58, we will use the Mies values until this matter has been resolved 
conclusively. We simply note here that reducing all of the (1,0) band Einstein 
A's by a common factor will reduce the g-factor by approximately this same 
amount. The Einstein A's are used to calculate oscillator strengths using 
eqtns. (12) and (14). The solar flux used for the (2,0) band is 4.7 X 10 13 
photons cm~^ s'VcnT* (Labs and Neckel 1968). Using these data, we have 
calculated g-factors for all the possible emission lines and these are shown in 
Figure 5a. The g-factor for the (1,0) band is given in Table 3. The value 
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* derived by Yamamoto (1982) is a factor of 1.8 larger than ours. Since we do not 
know the details of his calculation we cannot account for this difference. 

The fluorescence model of Schleicher and A'Hearn (1982) predicts directly 
the populations of the relevant rotational and vibrational levels of the ground 
state of OH. An individual line g-factor can be obtained from their results by 
multiplying the relative upper-state population by the Einstein A for that 
transition. The band g-factor derived from their data is a factor of ^2 smaller 
than our value (Table 3), and although it depends weakly on heliocentric 
velocity (through the UV pumping rate) this does not explain the discrepancy. 
Since we included only IR pumping in the (1,0) and (2,0) bands as excitation 
mechanisms, our calculation should yield a lower limit for the g-factors. 
Schleicher and A'Hearn (1982) neglected IR pumping mechanisms, and the IR band 
g-factors based on their model will give only the contribution made by UV 
excitation. By examining the relevant excitation and decay rates it is clear 
that direct IR pumping is about a factor of two stronger as an excitation source 
for IR emissions than is UV pumping. Schleicher and A'Hearn (private 
communication) have recently confirmed this by including IR pumping in both the 
(1,0) and (2,0) bands in their model. Their results for R = 1 a.u. and R ( = 
heliocentric radial velocity) = 0 km s"^ are shown in Fig. 5b. UV fluorescence 

causes population leakage from the ^k^/Z t0 ^1/2 ^ eve1s and > thererore * 
they predict the presence of twice as many lines as our calculation. Their new 
value for the (1,0) band g-factor is a factor of 1.5 larger than ours, of which 
2/3 represents direct IR pumping and 1/3 represents UV pumping. We note that 
the UV pumping of the IR (1,0) band represents an example of "Franck-Condon 
pumping", mentioned earlier (section II) as a potentially important excitation 
mechanism for diatomic molecules. Even in UH, which has a relatively narrow 
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* Franck-Condon envelope with relatively little vibrational pumping, the 
Franck-Condon pump contributes one third or' the total IR excitation. 

V. Discussion 

We have demonstrated that vibrational LTE cannot be maintained in the coma, 
except possibly within a very small region around the nucleus, and that 
scattering or' solar photons in fundamental ro-vibrational transitions of 
cometary parent molecules is expected to be the primary source of IR molecular 
emissions from comets. Moreover, these IR excitation rates are rather large 
(Table 3). The g-factor for the v-j band of COg exceeds the value of some of the 
strongest observed UV transitions (for example, the g-factor for the OH (0,0) 
band at x * 3090 % varies between 2 X 10‘ 4 and 10 X 10"^ photons s“ 5 ' molecule’ 1 -, 
depending on the heliocentric velocity). However, IR detectors are generally 
much less sensitive than UV or visible detectors (which can operate at the 
photon-counting limit), and large g-factors alone are no guarantee that the 
emissions are large enough to detect. The strengths of the IR molecular 
emissions from a typical moderately bright comet must be compared to the 
sensitivity limits of presently available IR instrumentation. 

We have chosen to predict the IR molecular emission from comet Bradfield 
(1979 X), a moderately-bright comet that was wel 1 -characteri zed by UV 
observations made with the International Ultraviolet Explorer satellite 
observatory in early 1980 (Feldman et al . 1980, Weaver et al . 1981). The 
observational parameters (R,a) are listed in Table 6, and are appropriate for 
comet Bradfield (1979 X) on 10 January 1980. The bright comet Bennett ( 1970 II) 
was observed by Ney (1974) in the IR (using filter photometry) under similar 


’ observational conditions during 1970. 
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Molecular production rates are given for comet Brad field (1979X) in 
Table 6. OH production rates were derived from measurements of its (0,0) band 
emission which was the brightest feature in the UV spectra (Weaver et al . 1981). 
Assuming that ^90% of the H^O molecules dissociate into OH molecules and H 


atoms, Weaver et al . (1981) found ii ,, = 2 X 1G 29 molecules s' 1 for v w n = 1 km 

, 'V n 2 

s . Since we use v^ ^ = 0.8 km s' 1 , we scale the production rate appropriately. 
The C0 2 production rate is obtained from measurements of the CO^ ion (Festou, 
Feldman, and Weaver 1981) and is approximately 10% of the H 2 ^ production rate. 

It is not yet clear which of COg or CO is more abundant in comets. Feldman and 
Brune (1976) measured a CO production rate which was ^0.3 of the H 2 U production 
rate for comet West (1976 VI). CO was tentatively identified in comet Bradfield 
(1979 X) (A'Hearn and Feldman 1980) and appears to be approximately five times 
less abundant than CO^ in this comet. We choose the CO production rate to be 
identical to that of C0 2 . The NH 3 production rate is the most uncertain 
quantity. Based on measurements of NH 2 in some comets, A'Hearn (1982) estimates 
the production rate of this molecule to be > 0.1% of the H,0 production rate. 
Since is thought to be a primary photodestruction product of NH 3 (Huebner 
and Carpenter 1979), NH 3 is expected to have a production rate similar to NH 2 * 
However, due to the uncertainties involved we have decided to use = 

nh 3 

0.1 H 2 0. Of course, if NH 3 is not in the nucleus then it will never be 
observed. 


■1 V, 


The average column density was calculated for each species within a 10" X 
15" field of view (FOV), the effective aperture for the IUE observations. These 
are listed in Table 6. An outflow velocity of 0.8 km s" 1 was used for all 
molecules except OH. The OH outflow velocity was assumed to be 1.15 km s _1 
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‘ (Weaver et al . 1981). The average surface brightnesses in the FOV were obtained 
using eqtn. (17) after scaling the g-f actors (Ta&le 3 and Figures 3, 5 and 6) to 
R * 0.71 a.u. using an inverse square law. We have listed the intensities in 
Table 6 using two systems of units: (1) Watts cm -2 sr" 1 , which is common in IR 

astronomy, and (2) Rayleighs, which is commonly used in UV and visible 
astronomy. The flux at the earth can be obtained from the intensity by 
multiplying by the solid angle of the FOV (3.5 X 10“ 9 sr, for our example). 

The predicted IR emissions are comparable to the brightest observed UV 
emissions. For example, the OH (0,0) band of comet Bradfield (1979 X) had a 
brightness of 370 kR, while the H Lya line brightness was 160 kR. Furthermore, 
we emphasize that the brightnesses of the IR emissions are not always 
distributed among many lines. There should be a tremendous reduction in the 
number of lines relative to the LTE case for H 2 Q, OH, and NH 3 emissions. Thus, 
high spectral resolution instruments may have an advantage in the detection of* 
these lines (we will elaborate on this later). 

The sensitivity of a presently available Fourier Transform Spectrometer 
(FTS), similar to the one used by Johnson, Fink, and Larson (1983) for their 
observations of comet West (1976 VI), is approximately 1.2 X IQ" 12 Watts cm" 2 
sr"* (this corresponds to ^210 Rayleighs near x = 2.7 ^m) in the region of the 
H 2 0 and OH bands (Larson, private communication). This sensitivity limit has 
been calculated assuming detector noise limited operation of the instrument, a 
spectral resolution of 4 cm”* (this is the same resolution used for the comet 
West (1976 VI) observations), a telescope aperture of 1.8 X 10 4 cm 2 , and an 
integration time of one hour. The FOV of their instrument is similar to -the 
value discussed above (depending on what detector they use, the FUV is either 
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- 20" X 20" or 10" X 10"), but for our calculation we assume their FOV Is exactly 
10" X 15". Comparing with Fig. 2, we see that at 4 cm -1 resolution individual 
lines of the band of H 2 Q are completely resolved, while there will be ^3-4 

lines of the V3 band of C0 2 in each resolution element (let us assume for the 

moment that the FTS is as sensitive near 4.3 u m as it is near 2.7 ^m, and that 
detector-noise limited operation is also possible at 4.3 iim). In Table 7 we 
give the average intensity in a single resolution element expected from the 
strongest HgO line in the v 3 band and from the sum of the three strongest C0 2 
lines in the p-branch of the C0 2 v 3 band. The expected signal to noise ratios 
are also given. These intensities appear to be easily detectable. 

The average optical depth in the FOV is estimated to be *Q.6 for the H^O 

line and ^0.4 for the C0 2 lines. While these average optical depths are not 
very large we must recognize that since the density distribution peaks very 
strongly towards the nucleus, radiative trapping and other optical depth effect 
may play a significant role in the inner coma. Optical depth effects are 
included approximately by calculating the average intensities using the model 
described in section III. These values are also listed in Table 7. While they 
are considerably smaller than the values derived assuming negligible optical 
depth, these intensities are still above the detector-noise limit of the FTS 
considered here. 

The specific intensity of the molecular lines can greatly exceed the grain 
emission. Using an empirical formula derived from the observational data on 
comet kohoutek (1973 XII) collected by Ney (1974), we have derived the average 
IR grain brightness in three wavelength regions for a "standard" comet at R = 
0.71 a.u. and A = 0.62 a.u. (see Table 8). The bandpass was chosen as 4 cm' 1 . 
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- The values measured for comet Bennett (1970 II) (R * 0.64 a.u., * * 0.76 a*u.) 
were scaled to a FOV of 10” X 15” using the scaling law suggested by Ney (1974) 
(Table 8). The average specific Intensity (Watts cnT 2 sr -1 /cnT l ) of the thermal 
grain radiation is compared with the line-center intensity of the strongest 
emission line for each molecule in Figure 7. If the IR molecular emissions are 
as large as predicted, then they should be easily seen in emission against the 
thermal continuum, using a sufficiently high spectral resolution instrument. 

Noise contributed by background thermal radiation from the sky, the 
telescope, or the instrument itself, may limit the detectability of these 
cometary lines by ground-based instruments. Since the background consists of 
greybody radiation at a temperature of ^270-300 K, it will be particularly 
troublesome in the vicinity of 10 urn. However, even for wavelengths as short as 
\ * 3 um, FTS instrument sensitivities may be limited by fluctuations in the sky 
brightness. Sky noise is difficult to predict theoretically and is highly 
dependent upon observing conditions (e.g., humidity, beam throw, etc.) and the 
particular instrument used. Background noise may be greatly reduced by severely 
restricting the bandpass of the instrument. For cryogenical ly-cooled 
space-borne instruments, photon noise from the cometary grain radiation may 
become important for low spectral resolution instruments (see Yamamoto, 1982). 

Telluric absorption is a particularly acute problem for ground-based 
observations because some of the molecules we are trying to detect in comets are 
also strong absorbers in the earth's atmosphere (e.g., HgO, C0 2 , and CO). The 
terrestrial absorption is not as severe for CO (Delbouille et al . 1981) and it 
may be possible to detect strong CO emissions from a comet whose emission is 
sufficiently Doppler-shifted away from the terrestrial absorption peak (the 
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* shift should probably be at least 0.3 cm -1 corresponding to 2 * 42 km s" 1 ). 
However, there appears to be little hope for ground-based detections of 
emissions from HgO, C0 2 » or the 0H( 1*0) band In the near IR. Detections of C0 2 
will probably only be accomplished from space, but it may be possible to detect 
H 2 o and OH IR emissions from comets using instruments at airplane altitudes. 

For example, the Kuiper Airborne Observatory (KAO) typically cruises at an 

altitude of 12 km, which is above 99% of the terrestrial wapor vapor. 

♦ 

Finally, we note that the high spectral resolution techniques (x/ax £ !0 6 ) 
of laser heterodyne spectroscopy might provide a relatively sensitive 
measurement of the NHj abundance in comets. We calculate that our instrument 
(Mumma et al . 1982) would be able to detect an NH 3 production rate of ^3 X 10 27 
molecules s" 1 for a comet at R * 0.71 a*u. and ^ » 0.62 a.u. (The heterodyne 
instrument has a diffraction-limited FOV which corresponds to *3" on a 60-inch 
telescope). The laser heterodyne receiver is also sensitive to the "hot" bands 
of CO g near 10 urn, as well as a host of emissions from other polyatomic 
molecules (although many of these have never been proposed as constituents of 
cometary nuclei). In fact, this method may be the only way in which a cometary 
parent molecule can be detected from the ground in the IR. 

VI. Conclusion 

It has been demonstrated that the physical conditions in the inner coma of 
a comet depart rather drastically from LTE. As a result, IR fluorescence of 
cometary parent molecules is expected to be a strong source of molecular line 
emissions. These emissions are predicted to be easily detectable using 
available instrumental sensitivities. However, due mainly to terrestrial 
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atmospheric absorption, the detection of IR molecular emissions will be a 
difficult experimental task and may require observations from airborne or 


space-borne platforms. Since detection of their IR molecular line emissions 
provides one of the few means by which cometary parent molecules can be 


identified directly (radio observations being the only other method which comes 
to mind), further research into this problem should be vigorously pursued. 


We acknowledge Dr. J. Bergstralh for deriving the empirical formula for the 
thermal grain radiation from a “standard" comet, 
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Comet 


Bennett (1970 II) 
West (1976 VI) 

Brad field (1979 X) 

a Al 1 values are derived 


Table 1 

Water Production In Comets 

Water Production Rate (molecules s" l ) a 


2 X 10 29 
4 X 10 29 
6 X 10 28 


from data presented In Festou et at. (1982) and are 


valid at R (a heliocentric distance) * 1 a.u. 
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Data for Collisions of H 2 0 with Various Molecules 


Molecule 


a Vf ( cm2 ) b 



co 2 

CO 

nh 3 


1.5 X 10“ 15 (v 2 ) 

3.5 X 10~ 16 ( v g ) 
2 X 10“ 20 (1,0) 

8 X 10~ 19 (v 2 ) 


2.5 X 10~ 15 

3.6 X IQ' 15 
3.2 X 10~ 15 
3.1 X 10“ 15 


Explanation of symbols and source of these data can be found in section III of 
the text. a + may be calculated from by invoking detailed balancing (see Sec. 
II, eqtn. 4). 

b The vibrational transition for the (T-V) process considered is listed in 


parentheses. 
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IR Molecular Emissions: g-Factors 


Molecule 

Transition 

Band Origin 

Band Strength 

Solar 

FI ux 

a Band 

9- 

•factor 



(cm* 1 ) 

(cm molecule* 1 ) 

(photons. cm" 2 
s'Ycnf l ) 

(photons s" 1 
molecule" 1 ) 

h 2 q 

v 2 

1595 

1.0 

X 

10- 17 

2.0 

X 10 13 

2.0 

X 

10" 4 

h 2 o 

v 3 

3756 

6.9 

X 

00 

1 

o 

3.9 

X 10 13 

2.7 

X 

9 

O 

r- 1 

C0 2 

v 2 

667 

8.3 

X 

10- 18 

8.6 

X 10 12 

7.1 

X 

IQ" 5 

C0 2 

v 3 

2349 

9.6 

X 

H 

o 

1 

H- 1 

2.8 

X 10 13 

2.7 

X 

10' 3 

CO 

(1*0) 

2143 

9.8 

X 

00 
r— a 
1 

o 

2.5 

X 10 13 

2.5 

X 

10~ 4 

NH 3 

v 2 

932,968 

2.2 

X 

o 

1 

*—• 

1.2 

X 10 13 

2.6 

X 

10 " 4 

OH 

(1*0) 

3570 

1.9 

X 

10" 18 

3.8 

X 10 13 

7.4 

X 

10" 5 


a Al 1 values are valid tor a heliocentric distance of 1 a.u. 
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Table 4 

Disequilibrium Distances 3 


Molecule: Transition 

n D (Debye) 

v (cm’ 1 ) 

A (S' 1 ) 

D f (km) 

(km) 

HgO; \>2 

1.88 

1595 

20 

11 

11 

v 3 


3756 

74 

- 

- 

rotation 


- 

4 X 10" 3 

- 

1000 

COgt vg 

0 

667 

1.4 

360 

20 

v 3 


2349 

399 

- 

- 

rotation 


- 

0 

- 

<9 

CO: (1,0) 

0.11 

2143 

34 

0.005 

0.03 

rotation 


- 

CD 

X 

o 

1 

GO 

- 

7.7 X 10 

NH 3 : vg 

1.47 

932,968 

15 

2.3 

0.3 

rotation 

j »» K, aJ=»1,aK*0 


m 

CM 

O 

»Hf 

X 

• 

220 

d » K, a-s 


- 

2 X 10’ 7 

- 

4.9 X ID 1 

J * K, aJ»l»aK 


«* 

0 

- 

o» 

a Di sequil ibrl urn distances are listed 

in the last 

two columns. 

HD is 

the 


permanent electric dipole moment. Other parameters are defined in section III 
of the text. 
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Table 5 

Photodestruct 1 on LI fetimes 

Molecule Lifetime ( T ,s) 



nh 3 

OH 


8.3 X 10 4 
5.0 X 10 5 

1.5 X 10 6 

5.6 X 10 3 

6.3 X 10 4 


a All values are valid at R (= heliocentric distance) ■ 1 a.u. Reference is 
Huebner and Carpenter (1979), except tor OH which is discussed In Weaver et al . 
(1981). 


46 


Table 6 




IR Molecular Emission from a Moderately Bright 
Comet at R ■ 0.71 a.u. and a * 0.62 a.u. 


Production 

Molecule Rate(mol s' 1 ) 


Column 3 Surface Brightness 3 

Density (cm' 2 ) Transition W cm" 2 sr' 1 Rayleighs 


h 2 o 

co 2 

CO 

nh 3 

OH 


1.8 X 10 29 
1.8 X 10 28 
1.8 X 10 28 
1.8 X IQ 28 
1.7 X 10 29 


3.2 X 10 15 

3.5 X 10 14 

3.6 X 10 14 

1.3 X IQ 14 
2.1 X 10 14 


v 3 

*3 

( 1 . 0 ) 

v 2 

( 1 . 0 ) 


1.0 X 10' 8 1.7 X 10 6 

7.0 X 10" 9 1.9 X 10 6 

6.0 X 10" 10 1.8 X 10 5 

1.0 X 10“ 10 6.7 X 10 4 

1.7 X 10' 10 3.1 X 10 4 


a Both the column density and the brightness are average values in a 10" X 15" 


aperture 
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Detection Capability for IR Molecular Emissions* 


Intensity (W cnT^ sr"*) 

Molecule Transition Optically thin S/N Optically thick 


S/N 


HgO v j band: 

°00 *01 


2.7 X 10' 9 2250 2.3 X I0* 10 


192 


COg v 2 band: 

P(12)+P(14)+P(16) 


9.5 X 10* 10 792 6.6 X 10* U 
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Assumptions are stated in the text (section V). Observational conditions are 
identical to those of comet Bradt'ield (1979 X) on 10 January 1980. Instrument 
sensitivities are those of the FTS described by H. P. Larson (private 
communication). The resolution was chosen to be 4 cm - ^. Detector-noise limited 
operation is assumed in the calculation of all signal -to-noi se ratios (S/N). 
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Table 8 



IR Background Radiation: Continuum Radiation from Cometary Grains 


Brightness (W cm” 2 sr” 1 ) 4 
3759 cm” 1 2347 cm” 1 1000 cm" 1 


“standard" comet at R-0.71 a.u. .. in Q 

and a- 0.62 a.u. 1.7 X 10” u 7.4 X 10” l ° 8.1 X 10”* 


comet Bennett (1970 II) 

at R»0.64 a.u. and o q « 

4-0.76 a.u. 1.7 X 10”* 8.0 X 10“* 7.0 X 10"° 


a These are average intensities in a 10" X 15" aperture. The values are 
integrals over a spectral interval of 4 cm" 1 . See section V of the text for 
further details. 
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Figure Captions 


Figure 1. The flux from comet Bennet (1970 II) , derived from the data of Ney 
(1974). The dotted curve has the shape of a blackbody spectrum at 
518 K, which Is the temperature most nearly describing the grains, 
according to Ney (1974). The excess cometary flux at *10 urn and *18 
urn are described as "silicate features" by Ney (1974). The co*tetary 
flux also shows an excess at the shorter wavelengths (x < 2.5 due 
to scattering of the solar flux by the grains. The solar flux shown 
Is the average continuum flux at 1 a.u. from the sun as derived from 
the data of Labs and Neckel (1968). The positions of the IR 
molecular emissions considered in this paper are indicated by 
arrows. 

Figure 2. Diagram showing the lowest energy levels of the (000) and (001) 

states of the H 2 0 molecule. The right hand side shows "ortho" levels 

which have a statistical weight three times that of the "para" levels 

located on the left hand side. The ortho and para levels are 

non-interacting. The energy levels are identified by J , where J 

**a 

is the total angular momentum quantum number and K a , K c are 
projections of J in the prolate and oblate symmetric top 
approximations, respectively. Other symmetry properties are 
indicated in the figure (see Herzberg (1945) for further 
explanation). The dashed lines show some of the rotational 
transitions in the ground state which will tend to relax the 
molecules to the lowest energy levels, as discussed in the text. For 
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completely relaxed HgO molecule? only the transitions marked with 
solid lines are Important* 

Figure 3. G-factors for Individual lines in the v 3 band of H 2 Q, the v 3 band of 
C0 2 » and the (1,0) band of CO. The values of the g-factors in 
absolute units (photons s" 1 molecule" 1 ’) are shown for a heliocentric 
distance of l a.u. Also given are the band g-factors, defined as the 
sum of g-factors over all lines in the band. 

Figure 4. Diagram showing the lowest energy levels of the ground vibrational 

state of NH 3 . Except for the K*Q levels, inversion doubling causes a 
splitting into symmetric (s) and antisymmetric (a) levels. This 
splitting (^ 0.8 cm" 1 ) is exaggerated in the figure for illustrative 
purposes. Total angular momentum quantum number, J, is given next to 
each level. K is the projection of J along the symmetry axis and its 
value is indicated along the horizontal axis. Positive and negative 
signs indicate the parity of the states. See Herzberg (1945) tor 
further details. Electric dipole allowed transitions which relax the 
NH 3 molecules to the J=K levels (as described in the text) are 
indicated in the figure. 

Figure 5. G-factors for individual lines in the v 2 band of NH 3 . The strongest 
transitions are identified by quantum number assignments. The values 
of the g-factors in absolute units (photons s" 1 molecule" 1 ) are shown 
for a heliocentric distance of 1 a.u. The band g-f actor is also 
given. The P, q, and R branches are shown separately and different 
wavenumber scales are used for the aQ and sq branches. 
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Figure 6. G-factors for the (1,0) band of OH. (a) shows our calculation which 
Includes only IR pumping, (b) gives the results of Schleicher and 
A 1 Hearn (private communication) including both IR and UV pumping. 
Individual line g-factors and their sum (the band g-factor) are shown 
for a heliocentric distance of 1 a.u. 

Figure 7. The intensity of the thermal radiation from the grains of a 

"standard" comet is indicated by the curve. The intensities at line 
center for the strongest lines in the bands considered in this paper 
are shown as "X". The heliocentric distance of tne comet is 0.71 
a.u., the geocentric distance is 0.62 a.u., and all intensities are 
average values in a 10" X 15" aperture. Our model predicts that at 
the position of the spectral lines, the flux collected from the IR 
molecular emissions will far exceed that of the grains in a 
sufficiently high resolution instrument. 
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SOLAR FLUX (photon* cm" 2 *“' / cm" 1 ) 
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